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Abstract 



This work investigates the existence and properties of a certain class of solutions of 
the Hill's equation truncated in the interval [r, r + L] - where L = Na, a is the period 
of the coefficients in Hill's equation, N is a positive integer and r is a real number. It 
is found that the truncated Hill's equation has two different types of solutions which 
vanish at the truncation boundaries r and r + L: There are always N — 1 solutions in 
each stability interval of Hill's equation, whose eigen value is dependent on the truncation 
length L but not on the truncation boundary r; There is always one and only one solution 
in each finite conditional instability interval of Hill's equation, whose eigen value might 
be dependent on the boundary location r but not on the truncation length L. 

The results obtained are applied to the physics problem on the electronic states in 
one dimensional crystals of finite length. It significantly improves many known results 
and also provides more new understandings on the physics in low dimensional systems. 
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1 Introduction 



Hill's equation is a class of second-order ordinary differential equations with coefficients 
of period a. The properties of solutions of Hill's equation has been investigated in detail 
by many mathematicians[f , 2]. It is well known that the eigenvalues and eigenfunctions 
of an ordinary second order differential equation depend on both the region in which the 
equation is given and the boundary conditions. If the equation is given in a finite interval 
and we are interested in a special class of the solutions - they vanish at the boundaries - we 
expect that both the eigenvalues and the eigenfunctions are dependent on the boundary 
locations and the properties of the finite interval. For an ordinary differential equation 
with periodic coefficients, in general it is usually expected that both the eigenvalues and 
the eigenfunctions are dependent on both the boundary locations and the interval length. 

By using the general theory of Hill's equation, this work investigates the existence and 
properties of this class of solutions of Hill's equation truncated in the interval [r, r + L] 
- where L = Na, a is the period of the coefficients in Hill's equation and iV is a positive 
integer and r is a real number. It is found that the truncated Hill's equation have two 
different types of solutions which vanish at the truncation boundaries r and t + L: There 
are always N — 1 solutions in each stability interval of Hill's equation, whose eigen value 
is dependent on the truncation length L but not on the truncation boundary r; There 
is always one and only one solution in each finite conditional instability interval, whose 
eigen value might be dependent on the boundary location r but not on the truncation 
length L. 

The Schrodinger differential equation on the electronic states in one dimensional 
crystals[3] can be considered as a special case of Hill's equation. Due to the fact that any 
real crystal always has a finite size, the electronic states in a real crystal of finite size can 
be better described as solutions of a special case of Hill's equation truncated in a finite 
length which vanish at the boundaries. The general results obtained here are applied to 
the physics problem on the electronic states in one dimensional crystals of finite length. 
It significantly improves many known results and also provides many new understandings 
on the physics in low dimensional systems. 
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2 Hill's equation 



We consider the Hill's equation having the following form: 

{p(x)y'(x)}' + {Xs(x) — q(x)}y(x) — 0, — oo < x < + oo (1) 

where A is a real parameter, p(x), q(x), and s(x) are real- valued and have the same 
period a. It is assumed that p(x) is continuous and nowhere zero and that p'(x), s(x) 
and q(x) are piecewise continuous with period a and there is a constant s > such that 
s(x) > s[2}. 

A periodic eigenvalue problem and a semi-periodic eigen value problem play a signif- 
icant role in the theory of Hill's equation (1)[2]: 

(i) . The periodic eigenvalue problem comprises (1), considered to hold in [0, a] and the 
periodic conditions: 

y(a)=y(0), y'(a)=y'(0). 

The eigenfunctions are denoted by ( n (x) and eigenvalues by v n (n = 0, 1,....). 

(ii) . The semi-periodic eigenvalue problem comprises (1), considered to hold in [0, a] and 
the semi-periodic conditions: 

y(a) = -y(0), y\a) = -y'(0). 

The eigenfunctions are denoted by £ n (x) and eigenvalues by /i n (n = 0, 1,....). 
The numbers u n and /i n occur in the order [2] 

< Mo < Mi < v x < v 2 < [ii < /x 3 < v 3 < is 4 < 

Let 0i (x, A) and (f>2(x,X) be the linearly independent solutions of (1), such that 
0i(O, A) = 1, 0i(O, A) = 0; 2 (O, A) = 0, 0' 2 (O, A) = 1, the real number D(X) defined by 
D(X) = 0i (a, A) +02(a, A) is called the discriminant of (1) and it plays a significant role 
in determining the properties of the solutions of (1). 

D(X) is an analytical function of A [2]. The values of A for which |-D(A)| < 2 form an 
open set on the real A-axis. This set can be expressed as the union of countable collection 
of disjoint open intervals (z^m, fi2m) and (/i2 m +i, ^2m+i)- Equation (1) is stable when A 
lies in these intervals thus the intervals are called the stability intervals of (1). The set 
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on the real A-axis with the stability intervals excluded can be expressed as the union 
of countable collection of disjoint intervals in which |-D(A)| > 2: (—00, u ], [//2m > A*2m+i] 
and [^2m+ij V2m+2\- Those intervals are called conditional instability intervals. Thus as A 
increases from —00 to +00, the conditional instability intervals and the stability intervals 
occur alternatively. 

3 The truncated Hill's equation 

A special class of solutions of the truncated Hill's equation are the main interest in this 
paper. We are trying to find out the eigenvalues A and eigenfunctions y(x), which are 
solutions of 

{p(x)y'(x)}' + {As(x) - q(x)}y(x) = 0, r < x < r + L (2) 

under the following boundary condition: 

y(r,A)=y(T + L,A) = 0, (3) 

where r is a real number, L = Na and N is a positive integer. 

Suppose y±(x, A) and y 2 (x, A) are two linearly independent solutions of (1), in general, 
the solution of (2) which satisfies (3) if it exists, can be expressed as 

y(x,A) = c 1 y 1 (x,A) + c 2 y 2 (x,A). r < x < t + L (4) 

In the following we try to find out the eigenvalues of (2) under the boundary condition 
(3) by assuming all non-trivial solutions of (1) are known. After each eigenvalue is found, 
the corresponding eigenfunction can be obtained easily. 

In principle, we should consider solutions of Equations (2) and (3) for A in (—00, +00). 
However, according to the Theorem 3.2.2 of [2], there is not a nontrivial solution of (2) 
satisfying (3) for A in (—00, z/ ]. Thus we need only to consider A in (z/ , +00). 

For A in (z/ ,+oo), D(\) may have five different cases depending on the value of A. 
Correspondingly, the two linearly independent solutions of (1) are also different: 
Case A. |£>(A)| < 2: 
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This is in the stability intervals of (1). Two linearly independent solutions of (1) can 
be expressed as [2] 

yi (x, A) = e ia ^ x Pl (x, A), y 2 (x, A) = e~ ia ^p 2 {x, A), 
where a (A) is a real number depending on A and 

< a(\)a < 7r 

and pi(x, A) and p 2 (x, A) have period a. All a(A) and p±(x, A) and p 2 (x, A) are functions 
of A. 

From (3) and (4) we have 

c 1 e^p 1 (r,\) + c 2 e~ taT p 2 (r,\) = 0, 

and 

Cl ^(r + L) pi ( T + LiX ) + c 2 e-^ T+L ^p 2 (T + L,X) = 0. 

Due to that 

pi(r, A) =pi(r + L, A), p 2 (r, A) = p 2 (r + L, A), 

we obtain either 

e ia(X)L _ e -ia(X)L = Q (A.l) 

or 

cipi(r, A) = anrf c 2 p 2 (r, A) = 0. (A.2) 

We can easily prove that neither pi(r, A) nor p 2 (r, A) can be zero. Suppose pi(r, A) = 0, 
we have pi(r + a, A) = and thus yi(r, A) = yi(r + a, A) = 0. Then according to the 
theorem 3.1.3 of [2], A must be in the range of [n 2m , H2m+i] or [^201+1,^201+2], in which 
I Z?(A) I > 2. This is in contradictory to the condition |-D(A)| < 2 here. Similarly p 2 (r, A) 
can not be zero. Thus we must have c\ = c 2 = and the non-trivial solutions obtained 
from (A.2) do not exist. 

Thus the existence of non-trivial solution of (2) and (3) requires 

a(\)L = j7r, j = 1, 2, 3, N-l. 
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Therefore for each stability interval, there are N — l values of A,-, where j = 1, 2, ....N — 1, 
for which 



Each eigen value for this case, is a function of L, the truncation length. But they all do 
not depend on the location of the truncation boundary r. 
Case B. D(X) = 2. 

There could be two possibilities [2]: 
Bl. (f) 2 (a, A) and (f)[(a, A) are both zero: 



This corresponds to v 2m+1 = ^ 2m + 2 . Two linearly independent solutions of (1) can be 
expressed as [2] 



and pi(x, A) and p 2 (x, A) have period a. It is always possible to choose C\ and c 2 , which 
are not both zero, to make y(r, A) = Ciyi(r, A) + c 2 y 2 (r, A) = and naturally we have 
y(r + L, A) = 0. Thus we have a solution A = A of (2) and (3) in this case and the 
corresponding eigenfunction is a periodic function. 
Obviously for this case, we always have 



a(Aj) = j ir/L. 



(5) 



</» 2 (a,A) = 0Ua,A) = 0. 



Vi(x, A) = pi(x, A), y 2 (x, A) = p 2 (x, A), 



y(r,X) 



y(r + a,A) 



0. 



B2. 2 (a, A) and <f)[(a, A) are not both zero. 

Two linearly independent solutions of (1) can be expressed as [2] 



yi(x,X) = px(x, A), y 2 {x, A) = A) +p 2 (^,A), 



and pi(x, A) and p 2 (x, A) have period a. 
Equations (3) and (4) become 



CiPi(r, A) + c 2 (rpi(r, A) +p 2 (t,A)) = 0, 



and 



c lPl {r + L,X) + c 2 ((r + L)p 1 (r + L,X) +p 2 (r + L,A)) = 0. 
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Thus, we must have, 

c 2 Lp 1 (r, A) = 0, 

and 

CiPi(t,A) + c 2 (tpi(t, A) + p 2 (r,A)) = 0. 

These two lead to 

Pl (r,A) = 0, and c 2 = 0. (52.1) 

Thus if we have a solution for this case, we must have (B2.1). 
Case C. D(X) > 2. 

Two linearly independent solutions of (1) can be expressed as [2] 

Vl (x, A) = e^ x Pl (x, A), y 2 (x, A) = e~^ x p 2 {x, A), 

where /5(A) is a real non-zero number depending on A and p\(x, A) and p 2 (x,\) have 
period a. 

Equations (3) and (4) lead to 

c 1 e^p 1 (r,\) + c 2 e- f3T p 2 (T,\) = 0, 

and 

c x e^ r+L) Vl {r + L, A) + c 2 e-^ T+L >p 2 (T + L, A) = 0. 

Due to 

Pi(r, A) = pi(r + L, A), p 2 (r, A) = p 2 {r + L, A), 

we must have 

CiPi(t, A) = 0, and c 2 p 2 (r, A) = 0. 

By the Sturm Separation Theorem[4], the zeros of pi(x, A) are separated from the zeros 
of p 2 (x, A). Thus pi(r, A) and p 2 (r, A) can not be zero simultaneously, either c\ or c 2 
must be zero. Whether we have a solution or not in this case is dependent on whether 
either one of 

Pl ( T ,X) = and c 2 = (C.l) 

and 

p 2 (r,A) = and ci = (C.2) 
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is true. From the discussion on (Bl), (B2), (C), we can see that in the conditional 
instability intervals in which D(\) > 2, if we have a solution of (2) and (3), we must 
have either z/ 2m+ i = V2m+2 or one of (B2.1), (C.l) and (C.2) when u 2m +i < v 2m+2 . Thus 
we always have y{r + a, A) = if y(r, A) = 0. Therefore for conditional instability 
intervals -D(A) > 2, the following equation is a necessary condition for having a solution 
of Equations (2) and (3): 

y(r + a, A) = y(r , A) = 0. (6) 

Case D. D(X) = -2. 

There could be two possibilities [2]: 

Dl. 2 (a, A) and (f>[(a, A) are both zero: 

2 (a, A) =<j)[{a,\) =0. 

This corresponds to /i 2m = A*2m+i- Two linearly independent solutions of (1) can be 
expressed as [2] 

yi(x, A) = si(x, A), y 2 (x, A) = s 2 (x, A), 

and Si(rr, A) and S2(^! A) have semi-period a: Si(x + a) = —Si(x). It is always possible 
to choose Ci and c 2 , they are not both zero, to make y(r, A) = Ciyi(r, A) + C2y2( r ) A) = 
and we naturally have y(r, A) = y(r + L, A) = 0. Thus we always have a solution of (2) 
and (3) in this case. The corresponding eigenfunction is a semiperiodic function. 
Obviously for this case, we have 

v(t,\) = y(r + a,X) = 0. 

D2. 2 (a, A) and 0'^a, A) are not both zero. 

Two linearly independent solutions of (1) can be expressed as[2] 

yi (x, A) = s ± (x, A), y 2 (x, A) = x s^x, A) + s 2 (x, A), 

and Si(x, A) and s 2 (x, A) have semi-period a. 
Equations (3) and (4) become 

cisi(r, A) + c 2 (tsi(t, A) + s 2 (r, A)) = 0, 
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and 

c lSl (r + L, A) + c 2 ((r + L)si(r + L, A) + s 2 (r + L, A)) = 0. 
Thus, we must have, 

c 2 Lsi(r, A) = 0, 

and 

CiSi(t, A) + c 2 (rsi(r, A) + s 2 (r, A)) = 0. 

These two lead to 

Si (t,A) = and c 2 = 0. (£2.1) 

Case E. £>(A) < -2: 

Two linearly independent solutions of (1) can be expressed as[2] 

Vl (x, A) = e^ x Sl (x, A), ^(x, A) = e~^ x s 2 {x, A), 

where /9(A) is a non-zero real number depending on A and si(x,X) and s 2 (x, A) are 
semiperiodic functions. 

Equations (3) and (4) lead to 

c 1 e /3 ^s 1 (r, A) + c^^^r, A) = 0, 

and 

Cie /3(A)(r+L) Si(r + L? A ) + c 2 e~^ T+L h 2 (r + L, A) = 0. 

Due to 

si(r, A) = (-l)%(r + L, A), s 2 (r, A) = (-1)%^ + L, A), 
we must have 

CiSi(r, A) = 0, anrf c 2 s 2 (r, A) = 0. 

By the Sturm Separation Theorem [4], the zeros of Si(x, A) are separated from the zeros 
of s 2 (x, A). Because Si(r, A) and s 2 (r, A) can not be zero simultaneously, either c\ or c 2 
must be zero. Whether we have a solution or not in this case is dependent on whether 
either one of 

Sl ( r ,A) = and c 2 = (E.l) 
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and 

s 2 (r,X) = and c x = (E.2) 

is true. From the discussion on (Dl), (D2), (E), it can be seen that in the conditional 
instability intervals in which D(X) < —2, if we have a solution of (2) and (3), we must 
have either \i 2m = ^2m+i or one of (D2.1), (E.l), (E.2) when a 2m < A^m+i- From these 
equations we can see that we also always have y(j + a, A) = if y(r, A) = 0. Thus for 
D(\) < —2, the same Equation (6) is a necessary condition of (4). 

It is easy to see that (6) is also a sufficient condition for having a solution of Equations 
(2) and (3) for |-D(A)| > 2: Equation (3) can be obtained by repeating (6) for N times. 
Thus the Equation (6) is a necessary and sufficient condition for having a solution of 
Equations (2) and (3) in a conditional instability interval for which \D(X)\ > 2. 

Equation (6) does not contain the confinement length L, thus the eigen value A in a 
conditional instability interval |-D(A)| > 2 might be dependent on r, but not on L. 

The Theorem 3.1.3 of [2] indicates that for any real number r, there is always one 
and only one A in each conditional instability interval [/x 2m , A^m+i] or [i^m+i, v 2m+2\, for 
which we have y(r,A) = y{r + a, A) = 0. We can label these A as A Tj2m and A T;2m+ i. 
According to the theorem 3.1.3 of e[2], the range of A Ti2m are in [u 2m , fi 2m+ i] and A r)2m +i 
are in [u 2m+1 , u 2m+2 ]. 

It is interesting to see how each of these r-dependent eigenvalues A changes as r 
changes. Here we only discuss A r 2m . A r 2m+1 can be very similarly discussed. 

If fJ>2m = n 2m+i, then we have the case Dl, the A r 2m will not change as r changes, 
we always have A r2m = ji 2m and the corresponding eigen function y(x,A Tj2m ) will be 
dependent on r, but it will always be a semi-periodic function. 

The more interesting case is when ji 2m 7^ At 2m+1 . According to the theorem 3.1.2 
of [2], i2m{x) and £ 2m+ i(:r) have exact 2m + 1 zeros in [0, a). Then according to the 
Sturm Comparison Theorem[4], the zeros of £ 2m (a;) and the zeros of £,2m+i( x ) must be 
distributed alternatively: There is always one and only one zero of £ 2m +i(:r) between two 
consecutive zeros of £ 2m (x), and there is always one and only one zero of £ 2m (a;) between 
two consecutive zeros of £ 2 m+i (••£)• 

We start from a zero Xi j2m of £ 2m (:r). If r = x^ 2m , we have a solution of (2) and (3) 
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K,2m = H2m, and y(x,A Tt2m ) = &m(x)- There is a nearby zero xi, 2m +i of &m+i(x). If 
t = xi,2m+i, we have a solution of (2) and (3) A Ti2m = /i 2m +i, and y(x, A T , 2m ) = 6m+i(^)- 
If we treat r as a variable and let r go continuously from a zero X\ j2m of £ 2m (:r) to a 
nearest neighbor zero £i, 2m +i of £ 2m+ i(:r), because A Tj2m is a continuous function of r[2], 
the corresponding solution of (2) and (3) A T)2m will also go continuously from /i 2m to 
/i 2m+ i. Similarly, if r goes from the zero a?i )2m+ i of £ 2m +i(:r) to the other nearest neighbor 
zero x 2i2m of £ 2m (x) continuously, the corresponding A r 2m will also go continuously from 
A^m+i to fi 2m . Since in the interval [0, a), both (, 2m (x) and (, 2m +i(x) have exactly 2m + 1 
zeros, in general A r 2m as function of r will always complete 2m + 1 ups and downs in an 
interval of length a. Similarly A Tt2m+ i as function of r will always complete 2m + 2 ups 
and downs in an interval of length a. 

Depending on the value of r, the solutions of (2) and (3) may have different forms. If 
r goes in one direction (either in the right or in the left) from a zero x lt2m of £ 2m (x) 
(y(x,A) = (,2m(x) here) to its nearest neighbor zero xi^m+i of £ 2m +i(x) (y(x,A) = 
^2m+i{x) now) then again to the next zero x 2j2rrt of £ 2m (a;) (y(x,A) = £ 2m (a;) again) con- 
tinuously, the corresponding A ri2m will also goes from /x 2m to fi 2m +i then again back to 
/i 2m continuously. But for r in the two open intervals (xi i2m , £i, 2m +i) and (xi^m+i, ^2,2m) 
of this path, the solution of (2) and (3) - the function y(x, A) (4) has two different forms: 
one has the form cie^^ 1 Si(x, A), and the other has the form c 2 e~^ K) x s 2 (x, A). In 
which one section it has which form is dependent on p(x), q(x) and s(x) in (1). 

A function with the form of Cie 13 ^ x si(x, A) or c 2 e~^ x s 2 (x, A), in which (3(A) > 0, 
is mainly distributed near either one of the two ends of the truncated region, due to the 
exponential factor. We can call it as an end solution of (2) and (3). End solutions may 
also exist in the instability intervals D(A) > 2 (f 2m +i < A < ^ 2m + 2 ) and have the form 
of either cie /3 ^ A ^ x pi(x,A) or c 2 e~^ A ^ x p 2 (x,A) in the truncated region, where pi(x, A) 
and p 2 (x, A) are period functions. These end solutions are introduced into the instability 
intervals when the truncated boundary r is not a zero of either one of the two eigen 
functions ( n (x) of the periodic eigenvalue problem or when r is not a zero of either one 
of the two eigenfunctions £ n (x) of the semi-periodic eigenvalue problem. That is the case 
|£>(A)|>2. 

In Figure 1 is shown A r as the function of r. In the open interval (^1,0,^1,1) and 
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the open interval (#1,1, £1,0 + a) the function y(x, A) (4) has different forms: one has the 
form cie 13 ^ x Si(x, A), and the other has the form c 2 e~ /3 ^ x s 2 (x, A), shown as a dashed 
line and a dotted line, indicating that the end solution is mainly distributed near either 
one of the two different ends of the truncated region. 

In summary thus there are two different types of solutions of the truncated Hill's 
equation: there are always N — 1 solutions in each stability interval, whose eigen values 
depend on the truncation length L but not on the truncation boundary r. There is 
always one and only one solution in each finite conditional instability interval \ji2m, A*2m+i] 
or [i/ 2m+ i, V2m+2\, whose eigen value might depend on the boundary r but not on the 
truncation length L. 

4 A special example: Electronic states in one dimen- 
sional crystals of finite length 

The one-dimensional Schrodinger differential equation[3] with periodic potential is a spe- 
cial case of Hill's equation (1), in which 

h 2 

P( x ) = 7^— > s ( x ) = !» Q( x ) = v ( x ) ( 7 ) 

is the periodic potential and A = E is the energy eigen value. Here k is the Planck's 
constant and m is the electron mass. It has been well known that the corresponding 
eigenvalues form energy bands and the eigenfunctions are Bloch waves [3]. The more 
general results obtained here for the truncated Hill's equation can be directly applied to 
the electronic states in one dimensional crystals of finite length. 

By using a Kronig-Penney model potential, Pedersen and Hemmer found that for a 
finite crystal of length L = Na, there are N — 1 states corresponding to each energy band 
n > 0, the energies of these A^ — 1 states map the energy bands [5] of the infinite crystal 
exactly and this mapping does not depend on the boundaries of the finite crystal. This 
is a special case of the more general result obtained in (5). 

By using the theorem 3.1.2 of [2], the author has given an analytical solution on 
the complete quantum confinement of one dimensional Bloch waves [6] in an inversion- 
symmetric potential. The results obtained in [6] greatly generalized and improved many 
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results obtained in [5]. 

The results obtained in this work can be directly applied to the electronic states in 
one dimensional crystals of finite length. The exact and general results are [7]: 
For one dimensional crystals bounded at r and r + L, there are two different types of 
electronic states: There are N — 1 states corresponding to each energy band of the Bloch 
wave, their eigenvalues A are given by (5), thus are dependent only on the crystal length 
L but not on the crystal boundary r and map the energy band exactly; There is always 
one and only one electronic state corresponding to each band gap of the Bloch wave, 
whose eigenvalue A might be dependent on the crystal boundary r but not on the crystal 
length L. Such a L-independent state can be either a constant energy confined band 
edge state (if r is a zero of a band edge wavefunction) or a surface state in the band gap 
(if r is not a zero of either band edge wavefunction). The results obtained in [6] are a 
special case of the more general results obtained here. 
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Figure Captions 

Fig. 1. A r as a function of r when /x 7^ A*i in interval Xi >0 < r < x\$ + a. The 
zeros of £o(x) are shown as solid circles and the zero of are shown as an open circle. 
However, in general, A r n as a function of r does not have to be linear. 
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